Topological Hysteresis in the Intermediate State of Type-I Superconductors 
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Magneto-optical imaging of thick stress-free lead samples reveals two distinct topologies of the 
intermediate state. Flux tubes are formed upon magnetic field penetration (closed topology) and 
laminar patterns appear upon flux exit (open topology) . Two-dimensional distributions of shielding 
currents were obtained by applying an efficient inversion scheme. Quantitative analysis of the 
magnetic induction distribution and correlation with magnetization measurements indicate that 
observed topological differences between the two phases are responsible for experimentally observable 
magnetic hysteresis. 

PACS numbers: 74.25.Ha,74.25.0p,89.75.Kd 



The structure of the intermediate state in type-I su- 
perconductors has a long history beginning with the pi- 
oneering papers of Landau 111 l3| a.nd continuing to the 
present day H 0, H H H OTlIS IH El Intermedi- 
ate state flux patterns closely resemble those found in a 
wide variety of hydrodynamic, chemical and solid state 
systems |^ J4J . Study of the intermediate state is 
therefore vital to a general understanding of pattern for- 
mation. Flux structures can be readily tuned with a 
magnetic field and imaged with magneto-optical (MO) 
techniques js], The observed patterns can then be 
correlated with underlying thermal, magnetic and resis- 
tive properties. Early MO images of the intermediate 
state revealed a variety of phenomena not predicted by 
the simple theory 0, 0] . The initial models were then re- 
fined to include domain branching and corrugation. Still, 
it is widely believed that a true thermodynamically sta- 
ble configuration of the intermediate state is the famous 
Landau laminar structure 0. 

Magnetic hysteresis is routinely observed in type-I su- 
perconductors and has generally been attributed to im- 
purities, grain boundaries, dislocations and other imper- 
fections of the crystal structure 0|. In this Letter we 
focus on the relationship between the topology of flux 
structures and their macroscopic magnetic properties. 
We find that a small residual hysteresis remains even 
in the most carefully prepared samples and present ev- 
idence that this hysteresis arises from the differences in 
the topology of the intermediate state between flux entry 
and flux exit. 

Samples were prepared from 99.9999% lead [i^ foils 
and rods. More than a dozen samples were prepared by 
using various annealing protocols or deliberately intro- 
ducing stress by cold rolling. The most reversible samples 
were obtained by melting lead between two Pyrex slides. 



Samples had thickness between d = 0.1 and 1.5 mm and 
were about 1.5 x 1.5 mm^ in planar dimension. The 
topological features described here were thickness inde- 
pendent above d ^ 0.5 mm, which indicates that they are 
not due to surface-related effects. We show data for sam- 
ples which had demagnetization factors of about N = 0.5 
(determined both from initial magnetisation and direct 
calculations [l^.) Quantum Design MPMS magnetome- 
ter was used for DC magnetization measurements. MO 
imaging was performed in a pumped flow-type optical 
''He cryostat using Faraday rotation of polarized light in 
Bi-doped iron-garnet films with in-plane magnetization 
15]. In all images the bright regions correspond to the 
normal state and dark regions to the superconducting 
state. 

Figurenshows typical magnetization loops for the cold 
- rolled sample (top panel) and the most reversible sam- 
ple of similar dimensions (bottom). The stressed sample 
shows considerably more hysteresis than the stress-free 
sample. The hysteresis increases for decreasing field in 
the stressed sample as expected from pinning, whereas 
the hysteresis disappears approaching _ff = in the re- 
versible sample, indicating complete Meissner expulsion. 
In that sample pinning is also absent at larger fields 
as well. Also shown in FigQ are MO images taken at 
the same temperature and magnetic fields indicated by 
arrows. The flux structure for the stressed sample is 
dendritic and a significant amount of flux is trapped at 
H = 0. In the stress- free sample the patterns are no- 
ticeably different, revealing flux tube phase upon flux 
penetration and well-defined laminar pattern upon flux 
exit. This behavior was observed at all accessible tem- 
peratures. 

Figure|21shows details of the evolution of flux patterns 
in the stress- free sample at T = 6 K. After cooling in zero 
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FIG. 1: Magnetization loops in stressed (top) and stress-free 
sample (bottom). Also shown MO images obtained at mag- 
netic fields indicated by arrows ( the numbers indicate images 
in succession). 



field, a full magnetization loop was measured with a max- 
imum magnetic field exceeding He ~ 260 Oe. Shown by 
the solid line, the M{H) loop exhibits magnetic hystere- 
sis at intermediate fields. The crucial question is whether 
this hysteresis is due to extrinsic factors (defects or resid- 
ual stress) or it is an intrinsic property of the interme- 
diate state. To clarify this, we performed zero-field and 
field - cooling experiments. The results are shown by the 
symbols in Fig. |5| The circles are obtained after cooling 
in zero field to 6 K and then increasing field, whereas 
squares indicate measurements after cooling the sample 
to 6 K in a particular field. If the hysteresis were due to 
pinning, the zero field cooled circles should coincide with 
the ascending branch of the M (H) loop but field-cooled 
squares should not. Instead, FigI21 clearly shows that 
both circles and squares coincide exactly with the directly 
measured magnetization loop, implying that the hystere- 
sis is not due to pinning. The fact that the ascending 
and descending branches merge at small and large fields 
is also inconsistent with pinning. Also, we observed no 
magnetic relaxation with either flux penetration or flux 
exit. These results strongly suggest that the hysteresis 
in stress-free samples is due to the topological difference 
between the closed flux tube phase and the open laminar 
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FIG. 2: Magnetization loop in a stress-free sample accom- 
panied by MO and zfc-fc measurements (see text). Circles 
are obtained after applying field after cooling in zero field. 
Squares show results of field-cooled measurements. 



phase. 

MO images shown in Fig. [5] reveal that after pure 
Meissner screening, the intermediate state appears, not 
as laminae, but as an assembly of normal tubes carrying 
magnetic flux and separated by superconducting regions. 
(MO image jjl shows the entire sample, others zoom in to 
reveal the structure) . These images show that flux tubes 
have a variety of structures - from simple monodomain 
to complex objects threaded with superconducting tubes. 
The flux tube phase favors hexagonal symmetry, almost 
exactly as modeled by Goren and Tinkham [l^. In all 
cases, this tubular phase has a closed topology that al- 
lows screening currents to circulate 0. Similar pat terns 
were directly observed in In 0, , Re , Sn [^l l22j and 
Hg iJj and it seems that closed topology tubular pattern 
is a generic feature of the intermediate state of pinning- 
free type-I superconductors upon flux penetration. 

Another possibility for the hysteresis is the edge bar- 
rier for flux penetration (inclu ding both, Bean-Livingston 
and geometric barriers) 0, 0|. However, such 
a barrier would result in delayed flux penetration and 
more negative values of magnetization compared to the 
thermodynamic values. Figure |21 shows that in our 
samples penetration occurs at the thermodynamic field 
Hc{l — N) Ri 120 Oe and magnetization at that point 
is as supposed to be in thermodynamic equilibrium, 
4ttM = —He, independent of the demagnetization factor. 
We attribute the weak influence of the edge barrier to the 
large thickness of our samples. On the other hand, the 
edge barrier could be involved in the formation of the 
observed topologies. When small fingers of the normal 
phase are formed at the sample edge, the surface barrier 
will prevent their continuous penetration into the interior 
and will break them into small flux tubes as suggested in 
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FIG. 3: T = 5 K, H = 100 Oe: (left) Distribution of 
the magnetic induction upon flux penetration (top) and exit 
(bottom); (right) corresponding patterns of shielding currents 
density obtained by numerical inversion. Intensity is propor- 
tional to the current density. (Color online) 




FIG. 4: r = 5 K, if = 300 Oe: (left) Distribution of the 
magnetic induction upon flux penetration (top) and exit (bot- 
tom);(right) corresponding contour plots of shielding cur- 
rents density obtained by the numerical inversion. Arrows 
show the direction of currents. (Color online) 



Ref.[T3. 

Furthermore, quantitative imaging of the magnetic in- 
duction can be used to visualize spatial distribution of 
shielding currents. This experimental information is im- 
portant for theoretical analysis involving current-loop 
models ^3 ^3 as well as for general understanding of 
pattern formation in type-I superconductors. We used 
recently developed fast inversion scheme [U]. Figure El 
shows the result obtained for H — 100 Oe. Left im- 
ages correspond to flux penetration (top) and flux exit 
(bottom) in our most reversible sample. The right panel 
shows corresponding distributions of the shielding cur- 
rents density (proportional to the brightness). Clearly, 
current distribution exhibits two topologically distinct 
patterns. The closed topology of small current loops vs. 
open topology of branched current streams. 

The direction of the currents flow is seen in Fig01 
which shows the reconstruction aX H = 300 Oe. At this 
field the features are larger and contour lines with di- 
rectional arrows can be used to better visualize the flow 
patterns. In both topologies, currents flow counterclock- 
wise - against the direction of Meissner currents flowing 
along the sample edges. At low fields Meissner currents 
dominate and hysteresis is negligible. It appears only at 
intermediate fields when the intermediate state consists 



of mobile ensemble of flux tubes on field entry and lam- 
inar structure that forms escape paths for flux expelled 
by the Meissner effect upon flux exit. 

To evaluate quantitative correspondence of the MO 
and M (H) measurements we calculate total mag- 
netic moment from the MO images by using AnAI = 
/ [H — B(r)]c?'^r. Magnetic induction is linearly propor- 
tional to the intensity, hence integrating images and using 
initial slope of the measured Ad{H) loop for calibration, 
M is obtained. Figure |5l shows that M{H) loop from the 
MO images (solid symbols) is in a good agreement with 
the direct measurements (open symbols). 

Observed topological hysteresis is also clearly seen on 
the proflles of the magnetic induction. Following Lan- 
dau j2| , it has generally been assumed that the magnetic 
field inside the normal phase in the intermediate state is 
close to mililla. Fi gure |S1 shows profiles of the 
magnetic induction measured at the same external field 
in the fiux tube phase (obtained on flux entry) and in 
a laminar phase (obtained on flux exit). While the field 
in the laminae is comparable to the critical field He as 
expected, the field above the flux tubes is much smaller. 
Indeed, the measurements are carried out w 10/^m above 
the sample, so the measured field is reduced compared to 
the values inside the tubes. However, simple numerical 
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FIG. 5: Comparison of DC magnetization measured by mag- 
netometer and reconstructed from MO images. Arrows show 
that at the minimum 4-kAI = —He as expected at thermody- 
namic equilibrium without edge barriers. 



FIG. 6: Profiles of the magnetic induction measured at T = 4 
K in the flux tube phase (flux penetration) and laminar phase 
(flux exit). Right panel shows corresponding MO images. 



analysis with appropriate dimensions does not reproduce 
such substantial reduction. It is possible that flux tubes 
widen when approaching the surface Q , but why it is not 
seen in the laminar phase? 

Alternatively, it is possible that when a flux tube ini- 
tially appears at the sample edge with critical magnetic 
field inside. Due to closed topology, the total magnetic 
flux in such tube is now conserved. The nucleated flux 
tube is driven by Meissner currents toward sample inte- 
rior (these currents flow everywhere on the surfaces per- 
pendicular to the magnetic held until it is stopped 
at the center or later by other tubes piling up from the 
center outward. When the flux tube reaches the inte- 
rior, its radius may increase to minimize the magnetic 
field energy. Real time imaging showed that flux tubes 
produced at the sample edge continue to travel toward 
the center and form an apparently outwardly expanding 
phase 22], as seen in the lower-right panel of FigEl The 
real-time observations were first made by Solomon and 
Harris in 1971 [l^. The rigorous evaluation of the free 
energy of even a single tube is not simple. The difficulty 
is the lack of a sharp interface between the superconduc- 
tor and tube interior, which decreases the surface energy. 
The tubes repel each other at the large distances due to 
interaction of screening currents, but they attract each 
other when their "cores" overlap. When tubes merge, av- 
erage magnetic field in the tube increases until it reaches 
He- At this stage, the honeycomb lattice is formed and it 
persists almost up to H — He- The observed topological 
hysteresis is observed only at the stage when magnetic 
field inside the tubes is less than He and before the for- 
mation of rigid hexagonal lattice, Fig|21 A detailed study 
of tube nucleation and expansion is needed to quantify 
this issue. 
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